A note on cluster methods for strongly correlated electron systems 
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We develop, clarify and test various aspects of cluster methods dynamical mean field methods using 
a soluble toy model as a benchmark. We find that the Cellular Dynamical Mean Field Theory (C- 
DMFT) converges very rapidly and compare its convergence properties with those of the Dynamical 
Cluster Approximation (DC A). We propose and test improved estimators for the lattice self energy 
within C-DMFT. 



The development of dynamical mean field methods has 
resulted in significant advances in our understanding of 
strongly correlated electron systems, in particular in the 
area of the Mott transition M. This method, captures 
the local effects of correlations such as the Kondo effect 
and the transfer of spectral weight between the coherent 
and the incoherent part of the spectral function. It suf- 
fers however from limitations arising from its single site 
mean field character such as the lack of k dependence of 
the self energy. In the context of disordered system, the 
dynamical mean field theory has a precursor in the fa- 
mous coherent potential approximation (CPA) for which 
cluster extensions such as the molecular CPA |7j have 
been formulated. Naturally, generalizations of these and 
other statistical mechanical approaches such as the Bcthc 
Peirls approximation to the area of quantum interact- 
ing systems, have been investigated recently M, ||, ||] 
i @ i H)j It ^ n ^ s area °f mves tig a tion is in its begin- 
ning stages, and comparative studies of the various meth- 
ods are important to increase our understanding of the 
strengths and the limitations of these methods, at a level 
comparable to our present understanding of the single 
site dynamical mean field theory. This note is a contri- 
bution in this direction. We focus, here on the CDMFT 
|| and the DCA Q method, because both methods have 
been proved to be manifestly causal, i.e. the output of an 
approximated solution of the cluster equations is causal, 
as long as a causal method is used for the solution of the 
impurity model. Because no numerical investigations of 
the C-DMFT have been carried out, we test the perfor- 
mance of this method in a simple soluble model that was 
introduced by I. Affleck and B. Marston S. It has a k 
dependent, albeit static, self energy, and therefore is a 
simple playground to explore the cluster method. 

In the first part of the paper we describe the CDMFT 
p|, and introduce a real space formulation of the DCA. 
to facilitate the comparison with CDMFT. In the second 
part we compare the predictions of DCA and CDMFT 
for the short distance behavior of correlation functions, 
in different cluster sizes against the exact solution. In 
the CDMFT approach the lattice self energy is a derived 
quantity which needs to be estimated from the cluster self 
energy, an auxiliary quantity which enters the dynamical 
mean field equations. In the second part of the paper we 



provide improved estimators for the lattice self energy 
and discuss how they improve the convergence to the 
exact answer as a function of the cluster size. 

Real space formulation the cluster schemes: A fairly 
general model of strongly correlated electrons contains 
hopping and interaction terms. It is defined on a lattice 
of L d sites in d dimensions, and we divide the lattice in 
(L/L c ) d cubic clusters of L d sites (more general forms 
can also be considered). We denote with e* the internal 
cluster position and with R„ the cluster position in the 
lattice (therefore the position of the i-th site of the n-th 
cluster is R„ + e, ) . The lattice Hamiltonian is expressed 
in terms of fermionic operators /^ a and /n m )y g and can 
be written as: 

H= Y. *«,/3( R »- R m)/]La/lW (1) 

n,a,m,fi,n' ,a' ,m' ,/3' 

a — i,a and a is an internal degree of freedom (i.e. a 
spin, spin-orbital or band index). Most cluster schemes, 
neglect the interaction terms between different clusters. 
The effects of those interactions, can be treated using the 
extended dynamical mean field approach [pj but we will 
not discuss these improvements in this paper. All the 
different cluster schemes can be formulated as a self con- 
sistent equation for the cluster self-energy which consist 
of the following loop:(i) start with a guess of the clus- 
ter self-energy (E c ) a ^, (ii) from the cluster self energy 
compute the Weiss function or host cluster propagator 
(Go)a,f3, which enters in the effective action for the clus- 
ter degrees of freedom, (iii) use the effective action com- 
pute the cluster Green function (G c ) a ,p, (iiii) compute 
the new cluster self-energy, (iiiii) iterate this loop un- 
til the convergence is reached. The DCA and CDMFT 
schemes differ in the way step (ii) is carried out. Within 
the CDMFT one obtains the Weiss function from the 
cluster self energy by the equation: 
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-^ (iu> n + n)l - t(K) - S c 



(2) 



where i(K) Q .^ is the Fourier transform of the hopping 
matrix in eq. d3) with respect to R„ — R m (k is a wave- 
vector in the Brillouin zone reduced by L c in each di- 
rection), bj n is the Matsubara frequency and /i is the 
chemical potential. Once the Weiss function has been 
computed one can obtain G c by functional integration of 
the single site action. Step (iiii) is carried out using the 
definition of the cluster self-energy: S c = G^ 1 — G^ 1 . 

To facilitate the comparison with the CDMFT, in the 
following we shall derive the DCA scheme using the real 
space formulation of the cluster. To lighten the nota- 
tion we will assume that the variable a is conserved to 
make all the cluster matrices diagonal in a and subse- 
quently we will drop this index. We take periodic bound- 
ary condition on the cluster and we define the matrix 
Ui.j(K) = exp(— iK ■ ei)5 ei . e . It is crucial for the follow- 
ing to note that the matrix i(K) has the following rep- 
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where k c are the cluster momenta. Therefore the matrix 
U (K.)i(K)W (K) is diagonal with respect to cluster mo- 
menta: (L/(K)t(K)L>T(K)) 4J = pi; kc e lkc ' (e '^ > t(K + 
k c ). Using this property one can write the (ii) DCA eq. 
in real space as: 
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(i)I-UtW(K) 



(3) 



Since the matrices UtW (K) and S c are diagonal with re- 
spect to k c , this equation coincides with the DCA equa- 
tions of Jarell et. al. after a Fourier transformation 
with cluster momenta. Once Go is known, G c is com- 
puted by functional integration of the cluster effective 
action and the new cluster self-energy is obtained by 
S c (ke j ) = ^ 1 (kc)-G- 1 (k c ). 

Eq. allows a direct formulation of DCA in real space 
and a detailed comparison with CDMFT. We also note 
this real space formulation can be used to defined many 
causal cluster schemes, by introducing a different matrix 
U(k) in the previous equation. 

A simplified one dimensional large-N model: compari- 
son between the exact solution and the predictions of the 
cluster schemes. In the following we focus on a simple 
one dimensional model, originally introduced and stud- 
ied by I. Affleck and B. Marston || in two dimensions. 
We compare the DCA and CDMFT schemes to its exact 
solution. This model is a generalization of the Hubbard- 
Heisenberg model where the SU(2) spins are replaced by 
a SU(N) spins, the on site repulsion is scaled as 1/N and 
the large N limit is taken. Its Hamiltonian reads: 



H - ~t/ Afi.o-fi+ha + fi+l,afi,<y) 



(4) 
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FIG. 1. x as a function of /3 for /j, 
The points are the exact solution 
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L c = 2, 5, 10 and xg DMFT for L c = 2, 3, 4. 
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where i = 1, . . . , L and a = 1, . . . , N and we take the 
large L and N limits. In the following we will use J as the 
unite of temperature and therefore we put J — I and we 
rescale the hopping term t — > t/J. The thermodynam- 
ics of this model can be solved exactly since in the large 
N limit the quantity x = j? E CT fi,^) fi+iA*) does not 
fluctuate. Indeed (||) reduces to a free-fermions Hamilto- 
nian with a "renormalized" hopping term t — ► t + \ and 
a self-consistent condition on %: 



X 



jY,f(PE k )cosk, E k 



-2(t + x) cosfc + /i 
(5) 



where /i is the chemical potential, f((3Ek) is the Fermi 
function and (3 is the inverse temperature. 
We now apply the DCA approximation to the Hamilto- 
nian (^) . As previously, the computation is simplified by 

the fact that the quantity X ci° A = jf E* //, CT (*)/i+i,<r(*). 
where i and i + 1 belong to the same cluster, does not 
fluctuate in the large-N limit. As a consequence the func- 
tional integral on the cluster degrees of freedom reduce 
to a simple Gaussian integral. Thus, imposing periodic 
boundary condition on the cluster, the eq. G~ 1 (k c ) = 
G~ 1 ( fc c) + S c(fc c ) implies: S c (fe c ) = 2 X ^ CA cosk c . Us- 
ing the second DCA equation which expresses the clus- 
ter Green-function as a function of E c (fe c ), we obtain the 
self-consistent DCA relation for x?i CA '- 



( f A = i^/P^Jcosfc, (6) 

-2t cos(k c + K)- 2Xci CA cos k c + fj, 



K,k c 
Ex,k c 



Note that in the infinite cluster limit one recovers the 
exact equation (g). 
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FIG. 2. Lattice self energy predicted by the different 
methods compared with the exact solution (triangles) for 
t-P-fj,— l. The dotted line is the result of the CDMFT 
estimator in [6], whereas the continuous line is the result of 
the first improvement discussed in this paper. The dashed, 
long dashed and dot dashed lines are respectively the DC A 
results for L c = 2, 5, 10. 

Now we focus on the CDMFT approximate solution. 
As in the DC A case, since the quantity (x^i DMFT )i = 
jf^,, fi a (t)fi+i :(J (t) does not fluctuate, one obtains: 



(£c)u = {xT Ubl )ihi-i + ixT M f T )ihi + i- This is 

the generalization of the corresponding DCA expression 
to a case without periodic boundary condition. Note that 
now the quantity {x^i DMFT )i may depend on the clus- 
ter index. Denoting the eigenvectors and the eigenval- 
ues of the matrix t + £ respectively ip"(K) and \ V (K) 
{v = 1, ...,L C ), the (ii) CDMFT equation, which ex- 
presses the cluster Green-function in terms of the cluster 
self-energy, reads: 



(G c 



L r 



^2^(K)(^(K)Y 
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Km 



+ H-\»(K) 



(7) 



Using this expression we finally get the self-consistent 
CDMFT equation on ixa DMFT )i- 



{xT MFT ) 



Lr 



^CW(CiW)7(»,), 



(8) 
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Notice that eq. H corresponds to the exact solution of a 
model defined by a Hamiltonian similar to (0) in which 
Ji y i+i equals 1 if i and i + 1 belong to the same cluster 
and zero otherwise. This implies in particular that in the 
infinite cluster limit the CDMFT approximation gives 
back the exact solution. We have numerically solved the 
self consistent equations (plflpP to compare the DCA and 
CDMFT predictions for different cluster sizes to the exact 
solution. 




FIG. 3. x as a function of /3 for /i = 1, and t = 1.. The 
points are the exact solution. The lines are, from top to bot- 
tom: x?i° A (dashed line), Xii CA (long-dashed line), Xc 
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(dot-dashed line), Xia (continuous line), \i 

(dotted line). Xi a DMFT1 IS the result obtained using the 
CDMFT estimator for the lattice self-energy in [6] whereas 
Xi a DMFT2 corresponds to the first improvement discussed in 
this paper. 

In Fig. (|l|) we plot the result of this analysis for 
£=l,/i = lasa function of j3 and for different clus- 
ter sizes. The two methods converge (the convergence 
is not uniform in f3) toward the exact solution for high 
enough L c but C-DMFT converges better that DCA. In- 
deed the C-DMFT result are already surprisingly good 
for L c = 2. However, as we shall discuss below, there 
are two different ways to compute x within the cluster 
methods. The one used here is based on a real space 
cluster intuition. The second one, relies on a momentum 
space intuition and computes the correlation functions 
from the k dependent lattice Greens function, this pro- 
cedure was advocated in ref and . we will show that 
it is also gives accurate results for small clusters. 

The lattice self-energy. We now address the compu- 
tation of the lattice self energy. In DCA a discretized 
form of the lattice self energy in momentum space, en- 
ters directly in the evaluation of Go- On the other hand, 
CDMFT focuses on estimating the cluster greens func- 
tion, and the lattice self energy does not participate in 
the mean field equations, and has to be estimated later 
from the cluster self energy. For the simplified large-N 
one dimensional model studied in this paper the DCA 
prediction for the lattice self-energy reads: 



yDCA 

n latt 



(k) = £(fc c ) = 2 X 



DCA COS k r 



(9) 



L a 



k c ,k c + j-]. Whereas 



L c 



where k belongs to [- 

for C-DMFT an estimator for the lattice self energy 
is constructed using the matrix SR n . a: i — <5_R„+ Q .i 
where a is the internal cluster index and i de- 
notes a lattice site. The simplest form m is: 



Z? a ? t MFT (k) = Ea,/j5t(fc)Sa 1 /9^(fc), where S a (k) is 
the Fourier transform of the matrix S with respect 
to the original lattice index i. S a (k) can be eas- 
ily written in terms of the matrix U(k) defined be- 
fore: S(k) a = U aya (k)/^J17 C . Therefore the relation- 
ship between the lattice and the cluster self energy 
reads: Zgf t MFT (k) = E a A W ( k )^U(k)]^/L c . For 
example, in the case of the two site cluster we find: 
Xg? t MFT (k) = X °i DMFT cos{k), whereas the exact solu- 
tion gives Yiiatt{k) = 2xexCos(k). As a consequence, even 
if the value of \ is wen predicted by the CDMFT there 
is a factor 2 between the two self-energies. The reason 
of this discrepancy may be understood writing the sim- 
ple estimator of the lattice self energy || in real space 
{^iatt)i-j = J2 a ,i3: a -i3=i-j(^c)a,p/L c . This means that 
the lattice self energy for a certain value of i — j is ob- 
tained averaging over all the cluster self energy elements 
corresponding to a — /3 = i — j. In the limit of an infi- 
nite cluster translation invariance implies that the clus- 
ter self energy coincides with the lattice self energy in the 
bulk. Therefore the factor 1/ L c cancels and we get the 
exact solution. However for a finite lattice there are only 
L c — 1 factors for i — j = 1, L c — 2 factors for i — j=2 , 
. . . L c — k factors for i—j — k. Therefore it is highly desir- 
able to have improved estimators for smaller size clusters 
in which the formula in which the average over all the 
factors having a — (3 — k is weighted by their number 
1/(L C — k). One could also think to put an extra weight 
to extract the lattice self-energy only from the sites in 
the bulk, for which the CDMFT result should be better. 
We propose new general class of estimators for the lattice 
self energy in terms of the cluster self energy, that inherit 
its causality property. 



V^latt) 



latt)i—j 



a,/3 : a—(3=i—j 



Wa,/3(Ec)a,f3 



(10) 



where the matrix w a _p is positive definite and 

J2 a ,f3:a-p=i-j w a,P ~* 1 for L c ~> °° ( tni s guarantees 
a good behavior in the infinite cluster limit). Using that 
the trace of the product of two positive definite matrices 
is positive, one can easily prove that if the cluster self- 
energy is causal this formula produce a lattice self-energy 
which is also causal. Note that (|l(j) does not change the 
behavior for an infinite cluster, but can really improve the 
results for finite cluster sizes. For example, in Fig. @ we 
compare the exact lattice self-energy to the DCA and the 
CDMFT predictions (using the initial estimator proposed 
in [pi and the simple improvement w a ^ = \j(L c — 1) 
which weights in the right way at least the terms with 
a — j3 = 1) for j3 = \i = t = 1. We remark that there 
is an excellent agreement between CDMFT and the ex- 
act solution after that our simple improvement has been 
taken into account already for L c = 2. Whereas the DCA 
prediction becomes good for L c > 5. 

Relation between lattice and cluster observables. Once 
the lattice self-energy has been obtained within a cluster 



method, the lattice Green function can be straightfor- 
wardly computed. This offers a different way of estimat- 
ing the first neighbor correlation function x, using the 
lattice Greens function. This quantity can be computed 
inside the cluster (xd) or using the lattice Green func- 
tion, obtained by the lattice self-energy, (xia) and the 
two results do not coincide in general. In the case of C- 
DMFT one can understand what are the approximations 
responsible for this difference and why they are small. 
The C-DMFT approach is based on the cavity procedure 
H which, if it was carried out exactly, it would give back 
the same answer for the lattice and cluster observables. 
However, in the approximated cavity procedure adopted 
in the C-DMFT, one assumes that the contribution to the 
effective action coming from tracing out all the degrees 
of freedom outside the cluster is purely gaussian. This is 
clearly not the case in general and it is the main reason 
for the non-coincidence of lattice and cluster observables. 

In Fig. (|) we compare the DCA and CDMFT predic- 
tions for the lattice and cluster values of \ to the exact 
solution for a two site cluster, for t = l,/x=lasa func- 
tion of (3. These curves display the typical behavior found 
also for other values of the control parameters: Xu^ A is 
quite better than its cluster counterpart, whereas the C- 
DMFT prediction is quite stable. This is probably the 
result of having an approximate cavity construction for 
the C-DMFT [§. Moreover we remark that X? a MDT2 i 
obtained using the first improvement for the self-energy 
discussed above, almost coincides with the exact solution. 
Comparing the C-DMFT and the DCA lattice values of 
X we note that C-DMFT gives usually a little bit better 
answer in terms of accuracy and convergence with respect 
to cluster size. This is perhaps due to the smoothness of 
the self-energy in the C-DMFT case. 

It would be nice to eliminate the cluster self energy 
altogether from the C-DMFT approach or to use a self 
energy without discontinuities in DCA, in the spirit of 
the work of Katsnelson and Lichtenstein Q. However, we 
were unable to prove manifest causality of this approach. 

In summary, in this short note we compared the per- 
formance of the DCA method with that of the cellular 
DMFT, in a very simple toy model. We have also pro- 
posed new estimators for the lattice self energy within 
CDMFT, which is more efficient. Our study shows that 
a direct application of CDMFT, i.e. without exploiting 
the flexibility inherent in the choice of basis in its most 
general formulation is very efficient in converging to the 
correct solution already for a two site cluster. Compar- 
ing the DCA and the C-DMFT predictions are a little bit 
better in terms of accuracy and convergence with respect 
to cluster size. DCA estimates of physical quantities, 
are most accurately carried out using the lattice Greens 
function, and not from the real space cluster correlation 
functions. This is stressed in ref H, who view DCA as a 
momentum space method. C-DMFT is not so sensitive 
to the choice of cluster or lattice estimators, because of 



the underlying a cavity construction of the derivation || . 
These results are very encouraging, and warrant further 
applications of these methods to more realistic and dif- 
ficult problems. Since the most glaring deficiency of the 
CDMFT method is that it does not attempt to take into 
account in a direct fashion the translation invariance of 
the problem, we concentrated on the phase of this model 
(finite temperatures) which is translationally invariant. 
We expect that CDMFT will perform even better for the 
ground state properties since in this case translation in- 
variance is broken by dimerization. 
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